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Abstract 

The multiplicity of inclusive photons has been measured on an event-by-event basis 
for 158- A GeV Pb induced reactions on Ni, Nb, and Pb targets. The systematics 
of the pseudorapidity densities at midrapidity (p ma x) an d the width of the pseudo- 
rapidity distributions have been studied for varying centralities for these collisions. 
A power law fit to the photon yield as a function of the number of participating 
nucleons gives a value of 1.13 ±0.03 for the exponent. The mean transverse momen- 
tum, (pt), of photons determined from the ratio of the measured electromagnetic 
transverse energy and photon multiplicity, remains almost constant with increasing 
Pmax- Results are compared with model predictions. 



1 Introduction 

The primary goal of ultra-relativistic heavy ion experiments is to study nuclear 
matter under extreme conditions, in which hadronic matter is expected to un- 
dergo a phase transition to a new state of matter, the Quark-Gluon-Plasma 
(QGP). For a thermalized system undergoing a phase transition, the variation 
of the temperature with entropy density is interesting as the temperature is 
expected to increase while below the transition, remain constant during the 
transition, and then increase again [1,2]. Temperature fluctuations have also 
been proposed as a signature of the existence of a tricritical point in QCD [?]. 
These behaviours can be studied by two experimentally measured quantities, 
viz., the mean transverse momentum, (pr), and the pseudorapidity density at 
midrapidity, p ma x, for varying impact parameter, or centrality, for a number of 
colliding systems. These variables also provide additional information to char- 
acterize the evolving system. p max provides a measure of the energy density 
which is important to understand the reaction dynamics [4,5]. In addition, the 
change in shape of the pseudorapidity distribution should be investigated in 
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detail because it may provide a clue to the formation of the QGP phase [6,7]. 

Except for a few measurements of photon [8-11] and neutral meson [12-15] 
distributions reported earlier, most studies have been restricted to charged 
particle measurements (a review of charged particle measurements can be 
found in [16,17]), due to the difficulty of precise measurements of photon dis- 
tributions. The inclusive photons provide a picture of the system at freezeout 
since the majority of the photons emitted from the reaction are decay prod- 
ucts of produced particles, viz., ttq and rj, (only a small fraction are emitted 
directly during the initial stage of the collision [18]). The shape and width 
of the pseudorapidity distributions of photons may be different compared to 
those of the charged particles. Photon multiplicity measurements have also 
become increasingly important because of the recent interest in simultane- 
ous measurements of the multiplicity of photons and charged particles in the 
search for Disoriented Chiral Condensates (DCC) [19,20]. The formation of 
a DCC is expected to give rise to large fluctuations in the relative number 
of emitted charged particles and photons, analogous to the Centauro and the 
anti-Centauro types of events observed in cosmic ray experiments [21]. Photon 
measurements can also be used to study flow [22] and intermittency behavior 
of events accompanying a possible phase transition. 

In this letter, we present the first measurement of the photon multiplicity and 
pseudorapidity distributions, together with the (p?) of photons produced in 
collisions of 158- A GeV Pb with Ni, Nb, and Pb targets, carried out in the 
WA98 experiment [23] at the CERN-SPS. 



2 Experimental Setup 

In the WA98 experiment, the main emphasis has been on high precision, si- 
multaneous detection of both hadrons and photons. The experimental setup 
consists of large acceptance hadron and photon spectrometers, detectors for 
charged particle and photon multiplicity measurements, and calorimeters for 
transverse and forward energy measurements. The present analysis makes 
use of the photon multiplicity detector (PMD), the midrapidity calorimeter 
(MIRAC) and the zero degree calorimeter (ZDC). 

The PMD was placed at a distance of 21.5 meters from the target. It was a 
preshower detector consisting of 3 radiation length (X ) thick lead converter 
plates in front of an array of square scintillator pads of four sizes, varying from 
15 mm x 15 mm to 25 mm x 25 mm, placed in 28 box modules. Each box 
module consisted of a matrix of 38 x 50 pads and was read out using one image 
intensifier + CCD camera system. The scintillation light was transmitted to 
the readout device by using a short wavelength shifting fiber spliced with a long 
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EMA (extra-mural absorber) coated clear fiber. The total light amplification 
of the readout system was ~40000. Digitization of the CCD pixel charge was 
done by a set of custom built fastbus modules employing an 8 bit 20 MHz 
Flash ADC system. Details of the design and characteristics of the PMD may 
be found in ref . [24] . The results presented here make use of the data from the 
central 22 box modules covering the pseudorapidity range of 2.9 < rj < 4.2. 

The MIRAC [25] was placed behind the PMD at 24.7 meters from the target. 
It consisted of 30 stacks, each divided vertically into 6 towers, each of size 
20 cm x 20 cm, and segmented longitudinally into electromagnetic (EM) and 
hadronic sections. The depth of an EM section is 15.6X (equivalent to 51% 
of an interaction length) which ensures essentially complete containment of 
the electromagnetic energy, with 97.4% and 91.0% containment calculated 
for 1 GeV and 30 GeV photons, respectively. Hadrons also deposit a sizable 
fraction of their energy in the EM section. The MIRAC measures both the 
transverse electromagnetic (EtjV 1 ) and hadronic (E!p d ) energies in the interval 
3.5 < rj < 5.5 with a resolution of 17.9%/ \/E and 46.1%/V^E, respectively, 
where E is expressed in GeV. The ZDC measures the total forward energy, 
E F) with a resolution of 80%/v'E' + 1.5%, where E is expressed in GeV. 



3 Data analysis 

3.1 Event Selection 

The data were taken during the December 1996 Pb beam period at the CERN 
SPS. The thicknesses of the three targets were 250 /xm, 254 /zm, and 213 (ttrn 
for Ni, Nb, and Pb, respectively. The fundamental "beam" trigger condition 
consisted of a signal in a gas Cerenkov start counter located 3.5 meters up- 
stream of the target and no coincident signal in a veto counter with a 3 mm 
circular hole located 2.7 meters upstream from the target. A beam trigger was 
considered to be a minimum-bias interaction if the transverse energy sum in 
the full MIRAC acceptance exceeded a low threshold. 

Beam pile up, where a second beam trigger occurs at a time when the detectors 
are integrating their signals from the triggered event, was rejected by (a) using 
the timing information from the trigger detectors, and (b) requiring that the 
sum of energies in the ZDC and the MIRAC were within 3a from the average. 
Downstream interactions were also rejected by requiring a coincident signal 
from the forward hemisphere of the Plastic Ball detector which surrounded 
the target. To correct for other sources of background, data were also taken 
with no target in place. The target-out contributions were negligible except 
for the most peripheral reactions. 
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The CCD readout of the PMD was cleared every 10 //s using a clear pulse of 

1 /is width generated every 10 /is. This ensured that there was no substantial 
noise buildup on the CCD pixels between successive event triggers. A gate of 

2 /is around the clear pulse was used to veto partially or fully cleared events. 
The clear clock operated asynchronously and was vetoed with a 5.6 ms wide 
pulse when a valid trigger occurred to allow complete readout of all pixels. 
A further check on possible pileup in the CCD cameras was made by using 
a 10 /is range TDC to measure the time difference between the arrival of the 
last clear clock and any valid event trigger. Events with multiple interactions 
within the 10 /xs between clear pulses were rejected in the offline analysis. 

The centrality of the interaction is determined by the total transverse energy 
measured in the MIRAC. For the (pr) analysis, which used the MIRAC data 
directly, the centrality was determined instead by the forward energy, Ep, 
measured in the ZDC. The centralities are expressed as fractions of the mini- 
mum bias cross section as a function of the measured total transverse energy 
or measured Ep. The most central selection corresponds to the top 5% of the 
minimum bias cross section and the peripheral selection corresponds to the 
lower 50—80% range. Extreme peripheral events in the 80—100% range were 
not analyzed. 

3.2 Data Reduction in the PMD 

The digitized pixel charges are processed by using a pixel-to-fiber map to 
form fiber signals corresponding to each scintillator pad. The signals from 
several neighbouring scintillator pads are combined to form clusters, charac- 
terized by the total ADC content and the hit positions. On average, there is 
a 92% probability for photons to shower in the lead converter and produce 
large signals. Compared to this, hadrons give a signal mostly corresponding to 
minimum ionizing particles (MIP). The majority of the hadrons are rejected 
by applying a suitable threshold on the cluster signal. A fraction of hadrons 
undergoing interaction in the lead converter produce signals larger than the 
threshold and appear as contaminants in the photon sample. The number of 
clusters remaining above the hadron rejection threshold is termed as '-/-like. 
The characteristics of the preshower PMD are described by the following two 
quantities [24]: 

= N^/N™ (1) 
f p = Nf^/N^ (2) 

where e 7 is the photon counting efficiency and f p is the fractional purity of 
the photon sample. iV* nc is the number of incident photons on the PMD and 
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N]i s is the number of photon clusters above the threshold. Both e 7 and f v are 
determined by a detailed Monte-Carlo simulation using the VENUS 4.12 [27] 
event generator with default parameter settings and the detector simulation 
package, GEANT 3.21 [28]. The details of the simulation can be found in 
Ref. [24]. No lower threshold on the energy spectrum of photons is applied 
in the simulation. The photon counting efficiency decreases with increasing 
hadron rejection threshold. The purity improves significantly with increasing 
threshold only up to ~ 3 MIPs and then rather slowly at higher thresholds. For 
practical purposes a 3 MIP threshold appears as an optimum choice for hadron 
rejection. With this value the photon counting efficiencies for the central and 
peripheral cases are 68% and 73%, respectively. The purity of the photon 
sample in the two cases are 65% and 54%, respectively. 

From the experimental data one determines N 1 -uk e , the number of clusters 
above the hadron rejection threshold. Using the estimated values of e 7 and f p , 
one obtains the number of photons incident on the detector in the event from 
the relation: 

iV 7 = Nj_ Hke ■ f p / e 7 (3) 



3. 3 Extraction of (pt) of photons 



In a given event, the average transverse momentum of produced photons may 
be expressed as 

M - % (4) 



where E^ 1 is the transverse component of the electromagnetic energy, and 
iV 7 is the number of photons in a given 77-region. In the WA98 experiment, 
E^ 71 and iV 7 are measured with the MIRAC and the PMD, respectively, on 
an event-by-event basis. These detectors have complete overlap in azimuth in 
the region 3.5 < rj < 4.0. Hence the data in this region are used for computing 
the (pt) using the above equation. 

In order to obtain the final E^ 1 for equation (4), the measured electromagnetic 
energy in the MIRAC towers must be corrected for (1) the hadronic contri- 
bution to the EM section of the MIRAC, and (2) the energy deposited in the 
lead converter of the PMD because of its position in front of the MIRAC. The 
final expression may be written as: 

ET = ElAET - fh ■ fjai ■ {E? ad /(l - f h )}} MlOi 

1 — f PMD 
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where, Ef m and E^ ad are the energies measured in the electromagnetic and 
hadronic sections of the MIRAC towers, fh is the fraction of the hadronic 
energy deposited in the EM section, f ba i is the balance factor taking into 
account the different responses for electromagnetic and hadronic particles in 
the EM section [25] , N is the number of towers in the MIRAC within the given 
i] range, 9i is the polar angle of the i th tower, and Jpmd is the fraction of the 
electromagnetic energy deposited in the lead converter of the PMD. The value 
of fpMD is found to be 15%. Details of the corrections are similar to those of 
Ref. [9]. 

3.4 Systematic Errors 

Several different sources contribute to the final systematic error in the determi- 
nation of the number of photons, N T The error due to the effect of clustering 
of the pad signals is the dominant one. This error is determined from the sim- 
ulation by comparing the number of known tracks on the PMD with the total 
number of clusters obtained after clustering. Apart from the effect of multi- 
plicity as discussed in [10], the arrangement of box modules in the present 
setup leads to splitting of clusters at the box boundaries. The net result is 
that the number of clusters exceeds the number of tracks with a deviation of 
3% in the case of peripheral events to 7% for high multiplicity central events. 

The number of 7-like clusters depends on the ADC value of the MIP peak 
as determined from Pb+Pb data. It has been estimated [24] that because of 
an admixture of ~20% photons in the MIP sample in the data, the extracted 
MIP ADC value is higher by 2 ADC channels. This causes the extracted 
photon multiplicity to be lower by 2.5%. We have included this as a source of 
systematic error. The error on e 7 because of the variation in pad-to-pad gains 
is found to be less than 1%. 

The purity factor, f p , depends on the ratio of the number of photons and 
charged particles within the PMD coverage. The systematic error associated 
with this ratio has been studied by using the FRITIOF [29] event generator 
in addition to VENUS. The average photon multiplicity by using FRITIOF is 
found to be higher by about 4% in peripheral and by 1% in central collisions 
compared to using VENUS. 

The combined systematic error on the final photon multiplicity is asymmetric 
and varies with centrality of the reaction. The total systematic errors are 
—3.4% and +7.5% for peripheral collisions and —7.1% and +3.4% for central 
collisions, varying little throughout the PMD acceptance. The negative error 
implies overestimation of number of photons. 

The photon counting efficiency determined in the present case relies on the 
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energy spectra of photons as given by the VENUS event generator. As the 
conversion probability for low energy photons falls sharply [26] with decreas- 
ing energy below 500 MeV, the estimate of e 7 may be affected if the energy 
spectra in the actual case is different. Preliminary measurements of the pho- 
ton energy spectra with the WA98 lead glass spectrometer indicate that there 
is an enhancement of photons below = 250 MeV/c over that given by 
VENUS. Taking into account this excess of low energy photons in the PMD 
acceptance, the photon counting efficiency would be overestimated by 2—9% 
for central events and 3—13% for peripheral events, the smaller value being 
for large pseudorapidity and the larger value being for the smaller pseudo- 
rapidity region of the PMD acceptance. The effect would be to increase the 
quoted PMD photon multiplicities. However, this uncertainty in the photon 
counting efficiency due to the uncertainty in the photon spectrum has not 
been included in the final errors. The WA98 low energy photon measurements 
will be described in a separate publication. 

The systematic error on the determination of (px) depends on the error in 
both E T m and iV 7 . The major source of errors in E T m are the contributions of 
hadronic energy deposited in the electromagnetic section of the MIRAC and 
the fraction of the electromagnetic energy deposited in the lead converter of 
the PMD. This is estimated to vary from 8.8% to 10.5% from peripheral to 
central events. Thus the combined systematic error on (pt) has been estimated 
to vary from 11.5% to 12.7%. 



4 Results and Discussions 

The photon multiplicity is determined using equation (3) on an event-by-event 
basis from the total number of 7-like clusters within the PMD coverage. The 
resulting minimum bias distributions of iV 7 are shown in Figure 1 for Pb+Ni, 
Pb+Nb, and Pb+Pb reactions at 158- A GeV. For comparison, the correspond- 
ing photon multiplicity distributions from the VENUS 4.12 event generator 
with default parameter settings, are superimposed in the same figure. The 
shape of these three distributions is governed by the collision geometry. For 
asymmetric collisions of Nb and Ni targets, small shoulders are present around 
iV 7 of 300 and 200, respectively. This shoulder is produced when a decrease 
in the impact parameter leads to little increase in particle production and the 
cross sections for these small impact parameters pile up at a fixed iV 7 . The 
VENUS event generator does not reproduce this shoulder well. The iV 7 values 
are higher for the data compared to those of VENUS. 

The pseudorapidity distribution of photons at different centralities are shown 
in Figure 2 (a) for Pb+Pb collisions. The data have been corrected for ge- 
ometry, efficiency, and purity factors. The filled symbols represent the mea- 
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sured data, and the open symbols are reflections of the filled symbols at 
Vc.m.{ = 2.95). The histograms show the corresponding distributions obtained 
from the VENUS event generator. The discrepancy between the VENUS re- 
sults and the data is about 10% for central collisions at midrapidity. The 
pseudorapidity distribution of photons at different centralities for Pb+Nb and 
Pb+Ni are shown in Figure 2 (b) and (c), respectively. The discrepancies be- 
tween the data and VENUS are larger for these reactions compared to that of 
Pb+Pb. 

The pseudorapidity distributions have been fitted with Gaussian distribu- 
tions. The peak position of the distribution (f] pe ak), the pseudorapidity density 
{pmax), and the width (a) are extracted from the fits. The i] pea k values for 
Pb+Pb distributions remain constant at 2.95 for all centralities. For Pb+Nb, 
Vpeak decreases from 3.03 ± 0.15 for central collisions to 2.95 ± 0.32 in the case 
of peripheral collisions. The corresponding values for the Pb+Ni reaction are 
3.10 ± 0.16 and 2.99 ± 0.31. Figure 3(a) and (b) show p max and o for the 
three reactions as functions of the number of participant nucleons, N part , at 
different centralities. The N part values are determined from the VENUS event 
generator. p max increases with N part while a doesn't change with increase of 

N pa rt ■ 

More insight into the systematics of the particle production can be obtained by 
computing the integrated number of photons (N* ot ) over the full phase space. 
This has been obtained from the Gaussian fit parameters to the pseudorapidity 
distributions. Figure 3(c) shows the extracted values of N* ot for Pb+Pb as a 
function of N part . The solid line shows a fit to the data using the function: 

N%* = C ■ (N part ) a (6) 

where C is a proportionality constant. The value of the exponent, a, is ex- 
tracted to be 1.13 + 0.03. To further explore the systematics, we have divided 
the full f] region (0—6) into two parts, one corresponding to the central rapid- 
ity region, 2.4 < rj < 3.4, and the other beyond this. For both of these cases, 
the Pb+Pb data yields a value of a = 1.13. In comparison, fitting the photon 
distribution from the VENUS event generator in the same two regions yield 
different exponents, with 1.10 ± 0.07 at mid-rapidity and only 1.0 ± 0.05 for 
the outer region. 

The mean transverse momentum, (px), as a function of p ma x, the pseudora- 
pidity density of photons at midrapidity, are shown in Figure 4 for Pb+Pb 
collisions. The data point at p max ~ 525 corresponds to the highest centrality 
bin, 0—1% of the minimum bias cross section for Pb+Pb. The systematic error 
on the absolute values are indicated by the upper and lower brackets on the 
data. The (pr) values are constant within the quoted error. For comparison, 
the results obtained from the VENUS event generator are superimposed in 
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the figure. The (pr) values obtained from VENUS are systematically higher 
compared to data, and show very little change with centrality. The indication 
of a small rise and saturation of (pr) seen in the data is similar to what has 
been reported for neutral pions [15]. 



5 Summary 

Photon multiplicities have been measured, on an event-by-event basis, in the 
forward region for 158- A GeV Pb induced reactions on Ni, Nb, and Pb tar- 
gets. The peak positions of the photon pseudorapidity distributions are found 
to shift forward in going from the Pb target to Nb and Ni, as expected. The 
photon pseudorapidity densities increase with the number of participant nucle- 
ons, while the widths of the pseudorapidity distributions remain constant with 
centrality. The integrated number of photons scales like N^ t1 almost inde- 
pendent of the rapidity range over which the integration is performed. This is 
similar to the predictions of the VENUS event generator, except that VENUS 
shows a smaller power for forward rapidities. The photon mean transverse 
momentum has been determined from the event-by-event ratio of the elec- 
tromagnetic transverse energy to the number of photons. After an initial rise 
with increasing p max , the value of (pr) is observed to remain constant. 



6 Acknowledgements 

We wish to express our gratitude to the CERN accelerator division for the 
excellent performance of the SPS accelerator complex. We acknowledge with 
appreciation the effort of all engineers, technicians, and support staff who have 
participated in the construction of this experiment. 

This work was supported jointly by the German BMBF and DFG, the U.S. 
DOE, the Swedish NFR, the Dutch Stichting FOM, the Stiftung fuer Deutsch- 
Polnische Zusammenarbeit, the Grant Agency of the Czech Republic under 
contract No. 202/95/0217, the Department of Atomic Energy, the Department 
of Science and Technology, the Council of Scientific and Industrial Research 
and the University Grants Commission of the Government of India, the Indo- 
FRG Exchange Programme, the PPE division of CERN, the Swiss National 
Fund, the INTAS under Contract INTAS-97-0158, ORISE, Research-in-Aid 
for Scientific Research (Specially Promoted Research & International Scientific 
Research) of the Ministry of Education, Science and Culture, the University 
of Tsukuba Special Research Projects, and the JSPS Research Fellowships for 
Young Scientists. ORNL is managed by Lockheed Martin Energy Research 
Corporation under contract DE-AC05-96OR22464 with the U.S. Department 



10 



of Energy. The MIT group has been supported by the US Dept. of Energy 
under the cooperative agreement DE-FC02-94ER40818. 

References 

[I] E.V. Shuryak and 0. Zhirov, Phys. Lett. B89 (1980) 253. 
[2] L. van Hove, Phys. Lett. B118 (1982) 138. 

[3] M. Stephanov, K. Rajagopal and E. Shuryak, Phys. Rev. Lett. 81 (1998) 4816. 
[4] J.D. Bjorken, Phys. Rev. D27 (1983) 140. 

[5] M. Kataja, P.V. Ruuskanen, L. McLerran and H. von Gersdorff, Phys. Rev. 
D34 (1986) 2755. 

[6] S. Sarkar, D.K. Srivastava and B. Sinha, Phys. Rev. C51 (1995) 318, Erratum- 
ibid. C51 (1995) 2845. 

[7] A. Dumitru, U. Katscher, J. A. Maruhn, H. Stocker, W. Greiner and D.H. 
Rischke, Z. Phys. A353 (1995) 187. 

[8] M.M. Aggarwal et al., WA93 Collaboration, Phys. Rev. C56 (1997) 1160. 

[9] M.M. Aggarwal et al., WA93 Collaboration, Phys. Lett. B404 (1997) 207. 

[10] M.M. Aggarwal et al., WA93 Collaboration, Phys. Rev. C 58 (1998) 1146. 

[II] R. Albrecht et al., WA80 Collaboration, Phys. Lett. B201 (1987) 390. 
[12] R. Albrecht et al., WA80 Collaboration, Nucl. Phys. A 488 (1988) 651. 
[13] R. Albrecht et al., WA80 Collaboration, Phys. Lett. B 361 (1995) 14. 
[14] R. Albrecht et al., WA80 Collaboration, Eur. Phys. J. C, 5 (1998) 255. 

[15] M. M. Aggarwal et al., WA98 Collaboration, Phys. Rev. Lett., 81 (1998) 4087. 

[16] J. Stachel and G.R. Young, Ann. Rev. Nucl. Part. Sc. 44 (1994) 537. 

[17] M.M. Aggarwal and S.I.A. Garpman, Int. J. Mod. Phys., E4 (1995) 477. 

[18] E.L. Feinberg, Nuovo Cimento A34 (1976) 391; 
E.V. Shuryak, Phys. Lett. B79 (1978) 135; 
B. Sinha, Phys. Lett. B197 (1987) 263; 

D.K. Srivastava and B. Sinha, Phys. Rev. Lett. 73 (1994) 2421. 

[19] J.D. Bjorken, K.L. Kowalski, C.C. Taylor, "Baked Alaska", SLAC-PUB-6109, 
1993 and ^ep-ph/9309235 ; 

K. Rajagopal and F. Wilczek, Nucl. Phys. B399 (1993) 395. 



11 



[20] M.M. Aggarwal et al., WA98 Collaboration, Phys. Lett. B420 (1998) 169; 
T.K. Nayak et al., WA98 Collaboration, Quark Matter '97 proceedings, Nucl. 
Phys. A638 (1998) 249c. 

[21] C.M.G. Lattes, Y. Fujimoto, and S. Hasegawa, Phys. Rep. 65 (1980) 151; 

Y. Takahasi et al., (JACEE Collab.), Proc. 7 th Int'l. Symp. on Very High Energy 
Cosmic Ray Interactions, (1992), Ann Arbor, Michigan, ed. L. Jones. 

[22] M.M. Aggarwal et al., WA93 Collaboration, Phys. Lett. B403 (1997) 390. 

[23] Proposal for a Large Acceptance Hadron and Photon Spectrometer, Preprint 
CERN/SPSLC 91-17, SPSLC/P260. 

[24] M.M. Aggarwal et al., A Preshower Photon Multiplicity Detector for the WA98 
Experiment, |hep-ex /9807026 , Nucl. Instr. and Methods (in press). 

[25] T.C. Awes et al, Nucl. Instr. Methods A279 (1989) 479. 

[26] M.M. Aggarwal et al., Nucl. Instr. and Meth. A372 (1996) 143. 

[27] K. Werner, Phys. Rep. C232 (1993) 87. 

[28] R. Brun et al., GEANT3 user's guide, CERN/DD/EE/84-1 (1984). 
[29] B. Nilsson-Almqvist, E. Stenlund, Comp. Phys. Comm. 43 (1987) 387. 



12 




Fig. 1. Minimum bias inclusive photon cross sections for Pb+Ni, Pb+Nb, and 
Pb+Pb reactions at 158-A GeV. Solid histograms are the corresponding distribu- 
tions obtained from the VENUS event generator. 
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Fig. 2. Pseudorapidity distributions of photons in Pb induced reactions at 158-^4 
GeV on (a) Pb, (b) Nb, and (c) Ni targets. The solid histograms are the correspond- 
ing distributions obtained from the VENUS event generator. 
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Fig. 3. (a) Pseudorapidity density (p m ax), (b) width of the pseudorapidity distribu- 
tions (a), and (c) integrated values of number of photons (iV* *), as functions of the 
number of participant nucleons at different centrality bins for Pb induced reactions 
on Ni, Nb, and Pb targets at 158- A GeV. The solid line in (c) is a power-law fit to 
the data, which yields the value of the exponent, a = 1.13 ± 0.03. 
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Fig. 4. The mean transverse momentum, (pr), of photons as a function of the 
pseudorapidity density of photons at midrapidity, p max , corresponding to different 
centralities. The {px) values obtained from the VENUS event generator for Pb+Pb 
are superimposed for comparison. 
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